TOWARD SIMPLIFICATION OF DYNAMIC SUBGRID-SCALE MODELS 

C. David Pruett* 

Department of Mathematics 
James Madison University 
Harrisonburg, VA, USA 22807 


Abstract 

We examine the relationship between the filter and the subgrid-scale (SGS) model for large-eddy 
simulations, in general, and for those with dynamic SGS models, in particular. Prom a review of the 
literature, it would appear that many practitioners of LES consider the link between the filter and the 
model more or less as a formality of little practical effect. In contrast, we will show that the filter and 
the model are intimately linked, that the Smagorinsky SGS model b appropriate only for filters of firafc- 
or second-order, and that the Smagorinsky model is inconsistent with spectral filters. Moreover! the 
Germano identity is shown to be both problematic and unnecessary for the development of dynamic 
SGS models. Its use obscures the following fundamental realization: For a suitably chosen filter, the 
compatible resolved turbulent stresses, properly scaled, closely approximate the SGS stresses. 


1 Introduction 

By definition, direct numerical simulation (DNS) is the numerical solution of the Navier-Stokes equations 
without recourse to empirical models. In concept, the fluid motions are resolved down to the Kolmogorov 
length scale, at which eddies succumb to viscous dissipation. In general, the computational workload for 
fully-resolved DNS scales as Re 3 , where Re is the Reynolds number. Consequently, for the complex, high- 
Reynolds-number flows of engineering interest, the computational requirements DNS are staggering and 
prohibitive. 

In contrast, in large-eddy simulation (LES), the larger scales of motion are resolved in space and time 
on a moderately coarse grid; however, the effect of the subgrid-scale (SGS) motions on the evolution of the 
larger scales is modeled. In practice, the decomposition into resolved and unresolved scales is accomplished 
by a spatial (temporal) filtering operation with an associated cutoff length (time) scale A. 

First introduced in the 1960’s, LES has experienced afxesurgence of interest since 1991, when dynamic 
SGS modeling was proposed by Germano and coworkers®. The advantages and difficulties associated with 
dynamic SGS models are now well established, and space does not permit elaboration. However, it is fair to 
say that the promise of dynamic modeling has not been fully realized largely because many of the proposed 
fixes to the shortcomings of dynamic models involve considerable additional complexity and computational 
overhead. 

Here, our purpose is to examine the connection between the choice of the filter and the subgrid-scale 
(SGS) model, with an eye toward the simplification of dynamic SGS models. 

From a review of the literature, it would appear that many practitioners of LES consider the link between 
the filter and the model more or less as a formality of little practical effect. Surprisingly little is written 
on this topic, Piomelli et al. 9 and Aldama 1 excepted. Specifically, regarding the conventional practice of 
LES, Piomelli et al. 9 observed: “In the past, however, the choices of model and filter have been regarded as 
completely independent.” Recognizing that the behavior of the SGS model strongly depends on the choice of 
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the filter, they attempted to address the issue of filter-model consistency on the basis of physical arguments 
and a priori tests, which involve comparisons of the exact and modeled SGS stresses computed by fully 
resolved DNS. Here, we approach filter-model consistency from a mathematical point of view. From the 
present analysis come a' number of revelations, some of which run counter tb conventional wisdom. 

In the next section, we discuss aspects of linear filters in general, focusing on their order properties. In 
the third section, starting from the Germanp identity and in the context of a dynamic SGS model, we derive 
a simple but accurate approximation for the. residual stresses. In the fourth section, we derive a similar result 
without first appealing to the Germano identity. In the fifth section, we briefly discuss the implications of 
these results. Conclusions are summarised in the final section. 


2 Linear Filter Operators 


Like differential operators, filter operators may be either continuous or discrete. Hare, for brevity, we 
consider only continuous filters (which are also referred to w *aoalog?); however, the con cl u si o ns drawn for 
contineus filters geaeraUae immediately to linear discrete fllters. Moreover, also for brevity, we consider only 
time-domain filters. However, the implications should apply to spatial filters as well. 




Let /(t,£) be a continous function of time and space, and let A, the “window” width, denote a charac- 
teristic time scale associated with the temporal linear filter Ff/(t), A]. As a specific example, consider the 


continuous, causal filter given by the integral equation 


F [/(*.*), A] 



/(r,f)dr 


( 1 ) 


FYom first principles of the Calculus, it is readily shown that lim&_M) F[/(t), A] = /{<). On the other hand, 
for a finite window A, the time-domain filter above tends to remove oscillations of high frequency relative 
to A while preserving low-frequency oscillations, which defines a “low-pass” filter. For applications to LES, 
we consider only low-pass filters. 

The effect of a filter is most apparent in Fourier space. Tb each filter is associated a transfer function 
H((l) that quantifies the amplitude and phase effects of the filter on oscillations of dimensionless frequency 
fl = <v4. For example, the transfer function associated with Eq. 1, shewn in Fig. I, is readily obtained 
by directly integrating F(e*“‘, A) for arbitrary w. Figure I reveals some undesirable traits of the filter: the 
anq>Htude decay is not monotonic, the amplitude envelope decays slowly (like I/O), and, consequently, the 
cutoff is gradual. 


Figure 1. Transfer function of integral filter of 
Eq. 1. 
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Again by analogy to difference operators, to each filter operator is associated an order-property that 
quantifies the behavior of the filter as 0 tends toward zero. The order property is revealed by the lending 
non-zero term in A in the Taylor-series expansion of the filter. For example, for a function of time only, the 
Thylor-series expansion with respect to A of the filter of Eqi Ife 

F[m, a] = m - 1 m +^rw + w 3 ) 



( 2 ) 



f 


A class of causal time-domain filters more suitable for LES than Eq. 1 is that of the so-called Butterworth 
filters. Figure 2 compares the moduli of the transfer functions of prototypical Butterworth analog (BA) low- 
pass filters of orders 1, 2, and 4, each of which has a nominal cutoff frequency ft' c = 1. The properties and 
design constraints of the first- and second-order low-pass BA prototypes can be found in Strum and Kirk 12 . 
The fourth-order BA prototype was developed by the author using Mathematica. The prototype BA filters 
shown in Fig. 2 are readily discretized and adapted to an arbitrary cutoff ft c . 


Figure 2. Comparison of Butterworth aqidog 
(BA) filters of orders 1, 2, and 4 with spectral 
(sharp cut-off) filter. 
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The order of a filter is closely related to the flatness of its transfer function at ft = 0. In general, for a 
filter of order n, £T(*)(ft)|a.o = 0 for all 0 <■ Jb < n, but H< n+1 )(0) ^ 0. In particular, Butterworth filters 
manifest several desirable properties for applications to LES: 1) stability, i.e., |ff(ft)| < 1; 2) their transfer 
functions decay monotonically with increasing ft, and 3) their transfer functions are maximally flat near 
ft = 0. For comparison, Fig. 2 also presents the idealized transfer function of an analog spectral filter, which 
can be considered of infinite order. 

For brevity in future discussions, let an overbar denote filtered quantities; that is, / 5 F[/(t), A] for 
some fixed A. From Figs. 1 and 2 it can be inferred that, except for the spectral filter, / / /. It is also 
clear that the spectral ideal is more closely approximated as the order n increases. As a result, high-order 
temporal filters are problematic for practical applications to LES, because they necessitate the storage of 
relatively more time history. This may be a principal reason that, to date, time-domain filtering has been 
avoided by practitioners of LES, as hinted by Moin and Jimenez 8 in their survey paper. 

One might naively assume (as did the author originally) that higher order is better. One of the more 
significant results of this work is to show that, in the context of LES, lower order filters are desirable for 
several reasons. This is particularly good news if one wants to consider the application of time-domain filters 
to LES, for example, as in Pruett 10 . 

To develop the results of the next sections, we make use of the Taylor series expansions of filter operators. 
For example, recall Eq. 1 above. More generally, assuming sufficient differentiability of tx, any time-domain 
filter can be expanded as 

u(t, x, A) = u + C\ Ait # + C 2 A^u ,; + c 3 AV" + ... (3) 

where primes denote temporal partial derivatives. For spatially multi-dimensional filter operators, similar 
expansions could be derived; however, their Taylor expansions would, of course, be multi-dimensional. In 
general, a filter is of order n in A provided Cn # 0, but c* = 0 for 1 < k < n. Thus, Eq. 1, for example, 
defines a first-order filter. Similarly, one can show that both Gaussian and tophat physical-space filters are 
first order. 


3 Conventional Dynamic SGS Modeling 

In tensor notation, the linearly filtered Navier-Stokes equations are given by 
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where repeated indices imply summation, and Ty is the SGS (residual) stress tensor defined as 

Ty = Uitfj - S?5j 


( 6 ) 


Equation 6 is exact; inexactness enters only when the Ty are modeled. We focus now on dynamic modeling 
of the residual stresses. 


Dynamic models 6 of the SGS stress tensor are rooted in theGermano® identity and typically exploit 
successive (spatial) “grid” and “test” filtering operations with associated length scales I and l, respectively. 
To transition to temporal filtering, we explicitly assume that j f f = r, where r > 0 is a parameter. 
Typically, r = 2. The Germano identity relates the resolved turbulent stress tensor £y and the “subgrid" 
and “subtest” stress tensors, Ty and Ty, respectively. Specifically, 


where 


and 


£y = lij Ty 


£ii = - fVfj 


(7) 


( 8 ) 


i! Ty 3 (9) 

The Germano identity is exact; moreover, its left-hand side is computible. It remains to model each of the 
terms on the right-hand Side, which is frequently accomplished via the Smagorinsky eddy- viscosity model, 
namely ? 

Ty ~ ^SijThk « 2Cl a |5|5y (10) 


Here, Jy is the Kronecker delta, 3y is the resolved-scale strain-ratafonsor, J5| m yfy , l is the char- 
acteristic length scale associated with the grid filter, and £ is the model constant of primary interest. In 
general, a criticism of eddy-viscosity models is their implicit assuiqptij&thaf the principal axes oftheresid- 
ual stress and the resolved-scale strain-fate tensors are aligned (MoW snd 'Jimenez*). Indeed, hi a recent 
experiment that examined turbulent three-dimensional boundar^hyer fi<w, Gompton and Eaton 4 found 
considerable misalignment between residual-stress and strain-rate tensors in the near*wall regfon. They 
concluded that eddy-viscbsiiy models are inappropriate for such flows. Moin and Jimmies 8 propose a more 
generally applicable model, for which each stress-tensor component has its own coefficient, namely 


Ty « (CikShj + Cy*?w)J a |3?| = r {CikPki + Cjkfiki) 


( 11 ) 


where 


.. ( 12 ) 

We will refer to Eq. 11 as the “generalised residual-stress model" The generalized residual-stress tensor is 
symmetric; hence, six independent coefficients mustbedetermined- Inprinciple, these coefficients can be 
uniquely determined by the dynamic procedure of Germano et al.®, as adapted to the generalized model by 
Moin and Jimenez 8 . We briefly describe the procedure below. 

It is now assumed that the subtest stresses can also be modeled by Eq. 11 with the same coefficients Cy , 
in which case 

Ty M (CikSkj + C jh Ski)P |3| (13) 

(In the author’s opinion, this assumption represents a considerable leap of faith, given that ry and Ty are 
not formal twins.) From Eqs. 7, 11, and 13, a set of integral equations for the coefficients Ca follows, namely 


where 


2£y » ( CikCtkj + CjkCtki - Cik$kj - Cjk0ki) 


ay 3 2P|5|5y 


(14) 

(15) 


We now exploit the Taylor-series expansion of the test filter to obtain a simple approximation to Eq. 14. 
For specificity, we use the filter of Eq. 1 and its expansion Eq. 2. Omitting details for brevity, we obtain the 
approximation 

2 Cn * (r s - 1) ( C ik 0 kj + Cjk&ki) + 2 L.O.E. + 0(l 2 A 2 ) = 2(r 3 - l)r 0 - + 2 L.O.E. + 0(1 S A 2 ) (16) 

where the leading-order error (L.O.E.) term in A is given by 


L.O.E. 


-7 [' 2 ( c “^ + c » i £ r ) - + ?£ $ r i ). 


We conclude that 


Ttj « 



L.O.E. 
r 2 — 1 


+ 0(1 2 A 2 ) 


(17) 


( 18 ) 


We now assume that the highest order term is insignificant. It remains to show that the second term 
(L.Q.E.) on the right-hand side of Eq. 18 is of lesser significance (on average) than the first. We assume 
here that all quantities have been previously scaled by appropriate reference values, so that we are dealing 
only with dimensionless quantities. It particular, lengths have been scaled by the wavelength of the largest 
eddies, and time has been scaled by the large-eddy turnover time. From Eqs. 8 and 2, to leading order in 
A, the resolved turbulent stresses are given by 


*-§($){£) ♦<**> /..<*»> 

In scaled variables, on the basis of reasonable assumptions and approximations (omitted for brevity), it dan be 
argued that the L.O.E. above is relatively unimportant whenever the dimensionless frequency 0 £ &&/<< 1 
(provided the filter is of first- or second-order, as will be shown). Specifically, for example, if rt c to 0.1, an 
entirely reasonable value in practice, then, Eq. 18 is approximated simply as 


T H 



( 20 ) 


Remarkably, Liu et al. 7 arrive at a result similar to Eq. 20 from experimented measurements in a turbulent 
jet. Specifically, they obtain several components of the SGS stress tensor of a jet by two-di m e n s io nal 
particle velocimetry. Whereas eddy-viscosity closures correlate poorly with the measured residual stresses, 
the resolved stresses correlate well. They propose the simple stress-similarity model 


Tij — c^Cij ( 21 ) 

where the coefficient ci is empirically derived. For r = 2, they obtain ci = 0.45 ± 0.15 for a “clipped” (no 
backscatter) SGS model by matching the exact and modeled SGS dissipation rates. For a model without 
clipping the optimal coefficient is approximately unity (Meneveau, personal communication). Either way, 
their result corroborates our observation that the residual and resolved turbulent stresses should be highly 
correlated. 

The implications of our present results for the practice of dynamic SGS modeling are both troubling 
and hopeful. Because the effects of numerator and denominator of the modeled residual stresses essentially 
“cancel” in the present analysis, dynamic SGS models, viewed in the present light, are ultimately independent 
of the form of their underlying model (whether Smagorinsky or otherwise)! This unanticipated result suggests 
that the whole concept of dynamic modeling needs re-examination. 

In hindsight, it appears to the author that the basis of dynamic models in the Germano identity is 
fundamentally flawed. The Germano “identity” is actually tautological, having been derived simply by 
regrouping and renaming certain quantities from the starting point This is not to say that the idea 

of dynamic modeling is flawed, only that there is no necessity for the Germano identity, as will be shown. 
Moreover, not only is the Germano identity unnecessary, it results in the practical difficulty associated with 
the vanishing denominator of the model coefficient. 



4 Alternate Approach to Dynamic Modeling 


In light of the discussion above, it is natural to ask: Can the residual stresses be modeled by the resolved 
turbulent stresses without appealing to the Germano identity 9 v 

By applying the general Taylor expansion Eq. 3 to Eq. 6, we obtain 

Tij = (c 2 - 2ca)uiuj A 2 + (C 1 C 2 - 3c 3 )(ui«" + t»Ju")A 3 + H.O.T. (22) 

where H.O.T. denotes higher order terms. Because the SGS-stress tensor ry arises solely from the quadratic 
nonlinearity of the NS equations, it is quadratic at leading order in A, provided that the filter is of either 
first- or second-order. On the other hand, if the filter is of order n > 2, then ry is of leading order n. Because 
the Smagorinsky model is of second-order in l (or equivalently, in A), it can be concluded immediately from 
Eqs. 10 and 22 that the model is appropriate only in the context of first- or second-order filters. Moreover, 
the use of Smagorinaky-based SGS models is totally inconsistent with spectral filters (which as we have said 
previously, can be considered of infinite order). This conclusion shopld hold regardless of whether filtering 
is accomplished in space or in the time domain. Our results are supported by experimental evidence. Liu 
«t al. 7 find High correlations between ry and Ay when filtering' is accomplished consistently with either 
Gaussian or physical-domain top-hat filters (both of which are first-order in our terminology). On the other 
hand, negligible correlations exist when a sharp cut-off filter is used in Fourier space (i.e., a spectral filter In 
our terminology). 

Of fundamental importance in dynamic models is the reared turbulent stress tensor £y (Eq. 8). As 
mentioned previously, the resolved turbulent stresses can hi 2 directly computed by filtering the resolved 
velocity fields u<. Let us now expand £y analogously to Eq, 22 above. To this end, we presume that the 
grid and test filters differ only in their respective filter widths. More precisely, if the grid filter is defined by 
Eq. 3, then the test filter is defined by 

u(t,x) = F £u(t,f), A] = u + ci(rA)ti' -I- C 2 (rA) 2 u" + c 3 (rA) 3 u"' + ... (23) 


with the same coefficients c< as in Eq. 3. From Eqs. 3, 8, and 23, and with the aid of Mathematics, it follows 
that 

£y * {<$ — 2o2)u<uJ(>*A)* +(cfr 2 — 2ctcjr 3 4- CiOjr 3 — + u^u”)A 3 4- H.O.T. (24) 

Rom comparisons of Eq. 22 and Eq. 24, we conclude that the SGS stresses can be approximated to leading 
order by the resolved stresses scaled by r 2 ; that is, 


Tij « 



How good is the approximation? Let jtbe approximation error E be defined 



(25) 

(26) 


from Eqs. 22 and 24, we obtain 

Eij = |3c 3 (r - 1) + *<*(3 - r) - cf] A 3 (« + utf) + H.O.T. (27) 

from Eqs. 22 and 27, we immediately conclude the following: ‘ ' 

1. If the filter is of either first- or second-order in A, then the approximation error is of higher order (3) 
than is the SGS stress (2), and the approximation is likely to be reasonably accurate (given additional 
constraints to be addrosed shortly). ; 

2. On the other hand, for any filter of order n > 2, the approximation error is of the same order as r 

itself; hence, r is likely to vanish in the noise of the approximation. 


( 28 ) 


3. If the filter is of order two (ci = 0.0, c<i ± 0.0) then 

Eij = 3c 3 (r - 1)A 3 (u'u" + u'u") + H.O.T. 

4. Ear from being inadmissible, as implied by the conventional dynamic modeling approach, r = 1 is 
optimal for second-order filters in that the leading error term vanishes. 

... 5. Because the residual stresses can be approximated directly from the resolved turbulent stresses, the 
Germano identity is unnecessary for the development of dynamic SGS models. 


5 Discussion 


Although, for brevity, the present results were derived using time-domain filtering, similar results could 
have been obtained for spatial filtering, albeit by more arduous mathematics. For example, Eq. 19 is the 
time-filtered analog to the space-filtered result of Clark et.al. 2 as interpreted by Speziale 11 , who reports that 



(!%l\ 

"' 3 12 \8xk) \dxk) 


+ 0( A 3 ) 


(29) 


Although Eqs. 20 and 25 are similar and in reasonable agreement for qioderately large r, they are pot 
identical. Whence the difference? Because the former originates from the Germano identity and the latter 
explicitly avoids it, we speculate that the discrepancy follows from the assumption that r^- and Ty of Eq. 7 
can both be modeled by formally identical models, despite some formal dissimilarity. 


If both analysis and experiment conclude that the residual stresses correlate closely with the (oom- 
putible) resolved stresses, then it is tempting to suggest for LES the use of SGS models that contain only 
scale-similarity terms. However, it is well known (e.g., Liu et al. 7 ), that scale-similarity models alone are in- 
sufficiently dissipative, and such calculations are almost guaranteed to blow up, particularly if the numerical 
scheme is non-dissipative. Our interpretation of the situation is as follows: the SGS models of LES must 
unfortunately play two roles: one physical and one mathematical. Whereas scale-similarity models appear 
sufficient to capture the physics of SGS energy transfer, additional dissipation (e.g., a Smagorinsky-like term) 
is necessary for mathematical reasons; i.e., to stabilize the numerical scheme whenever resolution is marginal. 
These roles are somewhat separated by mixed models (e.g., Bardina 2 ), which include both scale-similarity 
and dissipative terms. 


Although our results are completely consistent with the experimental results of Liu et a!. 7 , they are 
on|y partially consistent with the DNS results pf Piomelli et al. 9 , whose a priori tests show good agreement 
between modeled and exact stresses both for a mixed! model with a Gaussian filter and for the Smagorinsky 
model with a sharp cut-off filter. Whereas the former result is consistent with our findings, the latter is not. 
However, as Piomelli et al. 9 are careful to point out: "The fact that the SGS stress is essentially zero when 
the cutoff filter is used on the present [DNS] grid indicates that, with that filter, the grid may be capable of 
resolving the Reynolds stress and no model is needed.” Thus, the inconsistency may be more apparent than 
actual. We are currently conducting a priori tests to further validate our present analysis. 


6 Conclusions 

1. Mathematically tautological, the Germano "identity” is suspect as a basis for dynamic SGS modeling. 

2. A practical difficulty with dynamic SGS modeling, manifested in the vanishing denominator of the 
model coefficient, is directly attributable to the use of the Germano identity. 

3. The Germano identity is not only problematic, it is an unnecessary basis for dynamic SGS models. 

4. For first- or second-order filter operators, the computible resolved turbulent stresses, when properly 
scaled, closely approximate the residual stresses, without appeal to the Germano identity. 



5. In general, filters of higher than second order are in co n s i st ent with the Smagorinsky SGS model. 

6. In particular, spectral filters are inconsistent with the Smagorinsky SGS model. 

7. In LES, the SGS model plays two roles: one physical and one mathematical. To separate these roles, 
mixed models should be exploited. In mixed models, the scale-similarity term captures the physics 
and the dissipative term prevents numerical instability. .Common experience with LES reveals that the 
scale-similarity term alohe is insufficient. 

8. Tire scaling of the scale-similarity term of mixed models should depend on the choice of the parameter 
r relating grid and test filter widths. This has been overlooked in practice. 

9. A new model for the dissipative term, directly based on the computible resolved turbulent stresses, is 
sorely needed. 
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